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ABSTRACT 

The unprecedented sensitivity of Herschel coupled with the high resolution 
of the HIFI spectrometer permits studies of the intensity-velocity relationship 
I(v) in molecular outflows, over a higher excitation range than possible up to 
now. In the course of the CHESS Key Program, we have observed toward the 
bright bowshock region L1157-B1 the CO rotational transitions between J=5-4 
and J=16-15 with HIFI, and the J=l-0, 2-1 and 3-2 with the IRAM-30m and 
the CSO telescopes. We find that all the line profiles Ico(v) are well fit by a 
linear combination of three exponential laws oc exp(— \v /fo|) with vq = 12.5, 4.4 
and 2.5kms _1 . The first component dominates the CO emission at J > 13, 
as well as the high-excitation lines of SiO and H2O. The second component 
dominates for 3 < J up < 10 and the third one for J up < 2. We show that these 
exponentials are the signature of quasi-isothermal shocked gas components : the 
impact of the jet against the L1157-B1 bowshock (J\ ~ 210 K), the walls of the 
outflow cavity associated with Bl (Tk ~ 64 K) and the older cavity L1157-B2 
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(Tk ~ 23 K), respectively. Analysis of the CO line flux in the Large- Velocity 
Gradient approximation further shows that the emission arises from dense gas 
(n(H 2 ) > 10 5 - 10 6 cm" 3 ) close to LTE up to J=20. We find that the CO 
,7=2-1 intensity-velocity relation observed in various other molecular outflows is 
satisfactorily fit by similar exponential laws, which may hold an important clue 
to their entrainment process. 

Subject headings: ISM: jets and outflows — ISM: molecules — Stars: formation 

1. Introduction 

During the earliest protostellar stages of their evolution, young stars generate fast col- 
limated winds which impact against the parent cloud through shock fronts, generating slow 
"molecular outflows" of swept-up material. The intensity- velocity relationship Ico{ v ) ob- 
served in low- J CO lines in molecular outflows has been studied by various authors, as a 
possible test for discriminating between entrainment mechanisms. Downes & Cabrit (2003; 
hereafter DC03) showed that hydrodynamical simulations of jet-driven molecular outflows 
could successfully account for the observed relation Ico(v) in CO J=2-l. The sensitivity 
and the range of excitation conditions explored were somewhat limited, however. 

The heterodyne instrument, HIFI, onboard Herschefl now allows studies with unprece- 
dented sensitivity, of the dynamical evolution of gas in protostellar outflows and shocks at 
spectral and angular resolutions comparable to the largest ground-based single-dish tele- 
scopes (de Graauw et al. 2010). In particular, HIFI gives access to the CO ladder from 
J=5-4 up to J=16-15, probing a wide range of physical conditions. 

As part of the CHESS Key Program dedicated to chemical surveys of star forming 
regions (Ceccarelli et al. 2010), the outflow shock region L1157-B1 was investigated with 
Herschel. The protostellar outflow driven by the Class protostar L1157-mm (d= 250 pc; 
Looney et al. 2007) is the prototype of chemically rich bipolar outflows (see Bachiller et 
al. 2001 and references therein). Gueth et al. (1996) showed that the southern lobe of 
this molecular outflow consists of two cavities, likely created by the propagation of large 
bowshocks due to episodic events in a precessing, highly collimated jet. Located at the 
apex of the more recent cavity, the bright bowshock region Bl has been widely studied at 
millimeter and far-infrared wavelengths and has become a benchmark for magnetized shock 



1 Herschel is an ESA space observatory with science instruments provided by European- led principal 
Investigator consortia and with important participation from NASA. 
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models (see Gusdorf et al. 2008). Preliminary results (Codella et al. 2010) have confirmed 
the chemical richness of L1157-B1 and revealed the presence of multiple components with 
different excitation conditions coexisting in the Bl bowshock structure (Lefloch et al. 2010, 
Benedettini et al. 2012). 

In this Letter, we report on high-sensitivity CO observations with HIFI of L1157-B1, 
from J=5-4 up to 16-15, and complementary observations of the J=l-0, 2-1 and J=3-2 
with the IRAM 30m and the CSO telescope. 

2. Observations and data reduction 

2.1. The HIFI data 

CO transitions between J=5-4 and J= 16-15 were observed with HIFI at the position 
of L1157-B1 ajsooo = 20 h 39 m 10.2 s <5 J20 oo = +68°01'10.5". The observations were carried 
out in double beam switching mode. The receiver was tuned in double sideband and the 
Wide Band Spectrometer (WBS) was used, providing a spectral resolution of 1.1 MHz, which 
was subsequently degraded to reach a final velocity resolution of 0.5 km s -1 . The telescope 
parameters (main-beam efficiency r] m b, half power beamwidth HPBW) were adopted from 
Roelfsema et al. (2012; see Table 1). 

The data were processed with the ESA-supported package HIPE 6^| (Herschel Interactive 
Processing Environment). FITS files from level 2 data were then created and transformed 
into GILDAS§ format for baseline subtraction and subsequent data analysis. 

2.2. Complementary ground-based observations 

The CO J=3-2 line emission was mapped at the Nyquist spatial frequency across a 
region of 72" x 120" in the southern lobe of the L1157 outflow in June 2009 using the facility 
receivers and spectrometers of the Caltech Submillimeter Observatory (CSO) on Mauna 
Kea, Hawaii. Observations were carried out in position switching mode using a reference 
position 10' East from the nominal position of Bl. Small contamination from the cloud 
was observed, resulting in a narrow dip at the cloud velocity vi sr = +2.6 km s -1 (Bachiller & 



2 HIPE is a joint development by the Herschel Science Ground Segment Consortium, consisting of ESA, 
the NASA Herschel Science Center, and the HIFI, PACS and SPIRE consortia. 

3 http://www.iram.fr/IRAMFR/GILDAS 
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Perez-Gutierrez, 1997). The data were taken under good to average weather conditions, with 
system temperatures T sys in the range 750 — 850 K. An FFTS was used as a spectrometer, 
which provided a nominal resolution of 0.1 km s -1 . The final resolution was degraded to 
0.25 km s _1 . The final rms of the map is ~ 0.25 K per 0.25 km s _1 velocity interval. 

Deep integrations were performed at the frequency of the CO J=l-0 and J=2-l tran- 
sitions in June and August 2011 as part of an unbiased spectral survey of L1157-B1 at the 
IRAM 30m telescope (Lefloch et al. 2012, in prep). The EMIR receivers were connected to 
the 200 kHz resolution (FTS) spectrometers. Observations were carried out using a nutating 
secondary with a throw of 3', resulting in a narrow absorption feature at the cloud velocity. 

The observations (frequency, Obs_Id) and the telescope parameters are summarized 
in Table 1. Line intensities are expressed in units of antenna temperature corrected for 
atmospheric attenuation (for ground-based observations) Ta* ■ 

3. Results 

3.1. CO Spectral signatures 

The CO line profiles are displayed in Fig. 1 (J up > 5), Fig. 2a (J =1-0, 2-1) and Fig. 3 
(J=3-2) on a log-linear scale. This permits identification of three underlying components 
in the line profiles, denoted hereafter g±, g 2 and g^. Each component is well described by an 
exponential law Ico(v) = ico(0) exp(— |v/v |) showing the same slope at all J, but differing 
relative intensities. 

The high-excitation CO transitions J=13-12 and J=16-15 are well fit by the component 
gi alone, with v = 12.5 km s -1 (Fig. 1, top). In the lower J transitions, g\ still dominates 
the emission at high velocity v < — 20kms _1 (Figs. 1-2) and a simple scaling to the J=16- 
15 line profile in this velocity range allows us to determine the total g\ contribution to the 
integrated intensity of each CO line (Table 1). After removing the contribution of gi, the 
J=10-9 and J=9-8 line profiles appear to be well reproduced by the g 2 component alone, 
with vq = 4.4 km s" 1 (Fig. 1). Hence, the same procedure as above was applied to estimate 
the total contribution of g 2 across the CO ladder (Table 1). An emission excess with respect 
to the g 2 and g\ contributions is observed at velocities close to the cloud velocity at J < 7 
and actually dominates the emission in the low excitation transitions J < 2 (Table 1). This 
"residual emission" is well fit by the third exponential function g 3 oc exp(— |t>/2.5|), except 
for the transitions J < 2 for which v has a slightly steeper value 2.0. 

This profile decomposition is justified by the fact the the 12 CO line emission is optically 
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Fig. 1. — Montage of CO line profiles detected with H1FI. The line intensity is expressed 
on a logarithmic scale in corrected antenna temperature (T^*). The CO profiles are fit 
by a linear combination of two exponential functions g\ cx exp(— \v /12.5|) (blue) and gi oc 
exp(-|(w/4.4|) (red). 
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Fig. 2.— a) CO J=l-0 and J =2-1 line profiles observed toward L1157-B1 at IRAM. The line 
intensity is expressed on a logarithmic scale in corrected antenna temperature (T4*). The 
CO profiles are fit by a linear combination of three exponential functions g\ oc exp(— |u/12.5|) 
(blue) and g 2 oc exp(— \v /4.4|) (red) and #3 cx exp(— |t>/2.0)| (black). We superpose the fit 
to the full CO J =1-0 and ,7=2-1 emission with a dashed, magenta line, b) Similarity of the 
CO J=16-15 (E up = 750 K; thick black), H 2 3 12 - 3 03 (E up = 215.2 K; red dotted) and SiO 
J=8-7 line profiles (E up = 75.0 K; blue dashed), c) CO rotational diagrams at 20" spatial 
resolution for the three components : g\ (blue), g 2 (red), g 3 (black), d) x 2 distribution of 
LVG slab models for gi, g 2 and g^. Contours levels of 0.5, 1.0, and 2.0 are drawn with solid, 
dashed and dotted lines, respectively. 
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thin, as shown by comparison with 13 CO spectra, except at velocities very close to that of 
the ambient cloud for the low-J transitions. 

The fact that the slopes of gi, g 2 , and g 3 are independent of the CO transition considered 
is quite remarkable, and somewhat unexpected as one would naively assume the temperature 
gradients in shocked, accelerated gas to alter the shape of Ico(v) depending on the CO 
rotational level. Instead, it seems that the upper energy of the level only changes the relative 
importance of each exponential component in the resulting profile. In the following section, 
we show that this behavior is due to each component probing a distinct spatial region with 
almost uniform excitation conditions. 



3.2. Shock origin and physical conditions 

Figures 2a and 1 show that first g 3 , then g 2 and finally g\ dominate at progressively 
higher J up , which implies that the three components trace gas with progressively higher 
excitation conditions. We note that all three components emit over a wide range of veloc- 
ities and all display an emission peak at velocities close to the cloud velocity. Therefore, 
their excitation conditions cannot be determined from a simple analysis (e.g. line ratios) in 
different velocity intervals. Instead, we use our CO profile decomposition, which yields the 
total flux of each component as a function of J up (Table 1). The excitation conditions in 
each component are first obtained from a simple rotational diagram analysis of the HIFI and 
IRAM CO fluxes, after convolving to a common angular resolution of 20" (Fig. 2c). Further 
constraints on the kinetic temperature and density of the CO gas are then obtained using 
a radiative transfer code in the Large Velocity Gradient (LVG) approximation assuming a 
plane parallel geometry. We used the H 2 collisional rate coefficients of Yang et al. (2010) 
and built a grid of models with density between 10 4 cm -3 and 10 7 cm -3 and temperature 
between 10 K and 1000 K (250 K) to determine the region of minimum \ 2 as a function of 
density and temperature for g 1 (g 2 and g 3 ). We adopted a typical line width Av of 10 km s _1 
(gi), and 5kms _1 for (g 2 and g 3 ). 

3.2.1. The <7i component 

PACS observations of L1157-B1 have shown that the CO J=16-15 emission arises from 
a small 7 — 10") region, which peaks at « 5" North with respect to the nominal position 
of Bl, and is associated with a partly-dissociative J- type shock in the region where the 
protostellar jet impacts the cavity (Benedettini et al. 2012). We conclude that g 1 is the 
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spectral signature of this shock and refer to this gas component as the " g\ shock" in the 
subsequent discussion. 

As shown in Fig. 2b, the g\ component alone dominates the profiles of the H 2 3 i2 - 
3 03 (E up = 215 K) and SiO J=8-7 (E up = 75 K) lines, observed respectively with HIFI 
(Busquet et al., in prep) and the IRAM 30m telescope (Codella et al. in prep). In spite of 
large differences in their upper level energies, an excellent match is observed at all velocities 
between those tracers and the CO J=16-15 (E up = 751.8 K), except in the low-velocity 
range of the SiO J=8-7 transition, where a slight excess is observed. These three transitions 
are therefore probing the same physical region and the similarity of their line profiles justifies 
our use of the CO J=16-15 as a template for determining the contribution of g\ to each CO 
transition. 

Given the wide range of CO transitions and frequencies considered, the variations of 
the coupling of the telescope beam with the (off-centered) compact g± shock must be taken 
into account. This was done by convolving a fully sampled map of the SiO J=8-7 emission 
obtained at the IRAM 30m telescope (Codella et al., in prep) to the resolution of the HIFI 
beams of the CO transitions J=5-4 up to J— 8-7. The J=16-15 line flux in a 20" beam 
was directly obtained by convolving the PACS image of Benedettini et al. (2012) and the 
J=13-12 line flux was estimated under the assumption that the ratio of J=13-12/ J=16-15 
is preserved when degrading the resolution from 14" (10.7") to 20". This is supported by 
the fact that both lines show similar profiles. Since the IRAM J=2-l and HIFI J=16- 
15 observations have very similar angular resolution (Table 1), the J=2-l line flux was 
estimated under the same assumption as for J= 13-12. 

The rotational diagram of g\ is shown in Fig. 2c. The level populations in the HIFI 
range are well fit by a single rotational temperature T rot = 206 K and a beam-averaged 
column density iV(CO) = 2.3 x 10 15 cm~ 2 . The populations of the levels J up =l and J up =2 
lie a factor 3-7 above this trend, and could be the signature of a lower excitation component. 
LVG calculations were then carried out taking into account the CO line fluxes from J=5-4 
to J=20-19, using the HIFI and PACS data (see Benedettini et al. 2012). The results 
are shown in Fig. 2d in the form of \ 2 contours. The best-fit solution (x 2 = 0.35) is 
obtained for = 210 K and n(H 2 ) > a few 10 6 cm -3 , and a source-averaged column density 
iV(CO) = 0.9 x 10 16 cm" 2 for a source size of ~ 10". While our results are consistent with 
our previous PACS analysis (Benedettini et al. 2012), they favor dense solutions close to 
LTE with n(H 2 ) > 10 6 cm -3 and kinetic temperatures in the range 200 — 300 K (Fig. 2d). 
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3.2.2. The g 2 component 

Our CSO map shows that the CO J=3-2 outflow emission is dominated by the g 2 
component all over the Bl cavity,from the driving protostar L1157-mm down to the bowshock 
L1157-B1, at the cavity apex (Fig. 3). The exponent of g 2 remains constant at all the 
positions observed. This suggests that the g 2 component arises from the shocked gas in the 
walls of the Bl cavity. 

The excitation conditions in the g 2 component were derived by scaling the CO fluxes 
at each J to a common angular resolution of 20", following the same procedure as described 
above. The coupling between the telescope beam and the source was estimated from the 
convolution of our CO J=3-2 map to the resolution of the different HIFI beams. The data 
are well fit by a single rotational temperature T rot ~ 64 K and a beam-averaged column 
density iV(CO) = 4.0 x 10 16 cm" 2 . Our LVG calculations again favor an LTE solution with 
density above 10 5 cm -3 and kinetic temperature in the range 60 — 80 K (see Fig. 2d). The 
best-fit solution (x 2 = 0.27) is obtained for n(H 2 ) ~ 1.0 x 10 7 cm" 3 , T k = 64 K, and a 
source-averaged column density iV(CO) = 0.9 x 10 17 cm -2 for a typical source size of 20". 
Such a value of Tk is in good agreement with the value derived by Tafalla & Bachiller (1995) 
from multi-transition NH 3 observations using the VLA. The bulk of NH 3 emission arises 
from the cavity walls and peaks at the apex of Bl (Tafalla & Bachiller, 1995), supporting 
our interpretation of the g 2 component as tracing this cavity. 



3.2.3. The g 3 component 

The #3 component was identified as the residual emission in the CO line profile, in 
addition to the contributions of g\ and g 2 . The CSO data bring some insight into the spatial 
origin of this component. The CO J=3-2 emission from the older B2 cavity, south of Bl, 
is seen to follow the same intensity-velocity distribution Iqo^V) ^ exp(— \v /2.0|) as the g 3 
component observed toward Bl (see bottom panel in Fig. 3). We thus speculate that g^ is 
actually tracing shocked gas from the previous ejection, which led to the formation of the 
B2 outflow cavity. 

The excitation conditions in the g^ component were derived by scaling the CO fluxes 
at each J to a common angular resolution of 20", following the same procedure as described 
above. A rotational diagram analysis of the g% emission yields T rot = 26 K and a beam- 
averaged gas column density iV(CO) = 2.6 x 10 16 cm -2 . Our LVG calculations again favor 
an LTE solution with density > 10 5 cm" 3 and Tj< ~ 23 K (Fig. 2d). The best-fit solution 
was obtained for a source size of 25" and a source-averaged column density iV(CO) ~ 1.0 x 
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Fig. 3. — (Left) Southern outflow lobe of L1157 in CO J=l-0 (greyscale and black contours) 
and in SiO J=2-l (white contours) as observed at the PdBI (Gueth et al. 1996,1998). The 
HIFI band 1 FWHM (CO J=5-4) is shown as the black circle. Black squares mark the 
positions of the CO J=3-2 spectra plotted at right. (Right) Montage of CO J=3-2 spectra 
at various positions along the outflow between the protostar L1157-mm and the bowshock 
L1157-B1. The g% fit is drawn in red. The black solid line in the spectra of B2 and at offset 
position (0,-24) is representative of g^. The #3 component toward Bl (thick) is superposed 
on the spectrum of B2 to highlight their similarity. 
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10 17 cm -2 . The lower temperature of g 3 compared to g 2 is consistent with the B2 cavity 
being older than Bl, thus having experienced more post-shock cooling (Gueth et al. (1996) 
estimated an age of 3000 yr and 2000 yr for the outflow cavities associated with bowshocks 
B2 and Bl, respectively). 



3.3. The relation Iqo(v) revisited 




20 



|AV| (km/s) 

Fig. 4.— Plots of the observed I CO (v) in the J=2-l line for L1448, Orion A, NGC2071, 
L1551, and Mon R2 (open boxes and filled boxes; from Bachiller & Tafalla 1999). For each 
source, the fit Ico(v) oc exp(— |t>/t> |) is drawn as a black line, and the value of vo is given in 
km s -1 . 

Previous work described the relation Ico( v ) i n outflows by a broken power law, I co (v)oc 
v~ y with 7 ~ 1.8 up to line-of-sight velocities v break ^ 10 — SOkms" 1 and a steeper slope 
7 ~ 3-7 at higher velocities. This behavior is successfully reproduced by jet-driven flows, 
as a result of CO dissociation above shock speeds of 20kms _1 and of the temperature 
dependence of the line emissivity (see DC03). However, the underlying bowshock model 
predicts a power-law at low velocities on long time scales, rather than an exponential law, 
as observed in L1157-B1. It also predict a continuous range of temperatures in the swept-up 
gas, at odds with our finding that the L1157-B1 line profiles seem to be composed of three 
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quasi-isothermal spectral components. 

We show in Fig. 4 the CO J=2-l observations of five outflows previously studied by 
Bachiller & Tafalla (1999) and modelled by DC03 : L1448, Orion A, NGC2071, L1551, 
Mon R2. We display in dashed the best fit to the data with a single exponential Ico(v) 
oc exp(— \v/v \). A very good agreement is observed in all cases (Fig. 4), with values of v 
well in the range of those determined in L1157-B1. We conclude that an exponential relation 
Ico{ v ) ^ exp(— |f/fo|) is a good approximation to the observed intensity-relation not only 
in L1157-B1 but in molecular outflows in general, with a reduced number of free parameters 
compared to a broken power law. 

Our second main finding, that the exponential components in L1157-B1 appear quasi- 
isothermal and close to LTE, also has important implications. First, it shows that HIFI data 
are crucial to resolve ambiguities between sub-LTE vs LTE fits to CO excitation diagrams 
based on PACS data (Benedettini et al. 2012, Neufeld 2012). Second, it shows that the 
CO flux up to J up = 16 is dominated in each component by the densest, coolest postshock 
gas. Since such dense gas already reached a final constant speed, the broad velocity range of 
each exponential may require a broad range of view angles and/or shock speeds within the 
telescope beam. The exact origin of this spectral shape remains to be explained and may 
hold an important clue to entrainment and shock dynamics in molecular outflows. 
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Table 1. CO line observations. Telescope main-beam efficiency r] mh half-power beam 
width, rms (per lkms -1 velocity interval) and integrated fluxes of the shock components 
<7i, <?2 and g 3 are given. Temperatures are expressed in T A * units. 



Transition 


Frequency 




ObsID 


'Imb 


HPBW 


rms 


91 


92 


93 


Comment 




(GHz) 


(K) 






(") 


(mK) 


(Kkms- 1 ) 


(Kkms- 1 ) 


(Kkms" 1 ) 




1-0 


115.27120 


5.5 




0.78 


21.4 


2.0 


1.82 


13.8 


27.8 


IRAM 


2-1 


230.53800 


16.6 




0.59 


10.7 


2.7 


2.39 


29.7 


52.1 


IRAM 


3-2 


345.79599 


33.2 




0.65 


22.0 


130 




42.9 


17.5 


CSO 


5-4 


576.26793 


83.0 


1342181160 


0.75 


37.4 


8.0 


2.58 


39.8 


8.4 


HIFI 


6-5 


691.47308 


116.2 


1342207606 


0.75 


30.7 


5.0 


3.03 


45.3 


7.5 


HIFI 


7-6 


806.65180 


154.9 


1342201707 


0.75 


26.3 


7.1 


3.03 


36.9 


3.0 


HIFI 








1342207624 


0.75 


26.3 










HIFI 


8-7 


921.79970 


199.1 


1342201554 


0.74 


23.0 


10.0 


4.55 


27.4 




HIFI 








1342207323 


0.74 


23.0 










HIFI 


9-8 


1036.91239 


248.9 


1342200962 


0.74 


20.5 


7.7 


4.09 


23.8 




HIFI 








1342207641 


0.74 


20.5 










HIFI 


10-9 


1151.98544 


304.2 


1342207691 


0.64 


18.4 


36 


3.79 


9.61 




HIFI 








1342196511 


0.64 


18.4 










HIFI 


13-12 


1496.92291 


503.2 


1342214390 


0.72 


14.1 


46 


4.55 






HIFI 


16-15 


1841.34551 


751.8 


1342196586 


0.70 


11.5 


26 


2.27 






HIFI 



